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Spin liquid has been a cornerstone in the gauge theory 
approach for strongly correlated electrons^ as Fermi liq- 
uid in the Landau-Ginzburg- Wilson framework for phase 
transitions 2 -. An effective field theory is often given by 
compact U(l) gauge theory-^, implying that "magnetic 
monopole" excitations should be irrelevant in order to be 
self-consistent for the theory. Two kinds of mechanisms 
have been proposed, resulting from either spinon dynam- 
ics with Fermi surface^— or existence of a topological 
term£ associated with anomaly in the Dirac theory 9, 1°, 
After deconfinement is demonstrated, an important task 
is to solve the non-compact U(l) gauge theory. 

A standard technique is the large-A approximation, 
where the spin degeneracy of a spinon is extended from 
(7 =t,| to a = 1, ...,N. The N — > oo limit was believed 
to suppress higher-order quantum loop-corrections in the 
Fermi surface problem-^ just as the case of the relativis- 
tic invariant theory^. Recently, it was clearly demon- 
strated that the Fermi surface problem is still strongly 
interacting even in the large-iV limit, meaning that all 
planar Feynmann diagrams should be summed as the 
non-abelian gauge theory with Lorentz invariance 12 . 

This observation suggests that dynamics of fermions 
(spinons) and gauge fluctuations can be modified by 
more-loops quantum corrections, that is, the exponent in 
the frequency dependence of the spinon self-energy may 
have a nontrivial correction, affecting transport proper- 
ties of spinons. Actually, the lowest-order vertex cor- 
rection associated with the Aslamasov-Larkin diagram is 
shown to cause such a correction proportional to 1 /N al- 
though the lowest-order vertex correction associated with 
the ladder diagram does not change the dynamics of both 
fermions and collective bosons^. 

In this paper we perform an infinite-order summation 
for the ladder-type vertex correction and find no anoma- 
lous correction for the exponent in the frequency depen- 
dence of the fermion self-energy. In other words, dynam- 
ics of both spinons and gauge bosons remains the same as 
the case without vertex corrections^. We prove this re- 
sult based on the Ward identity-^, asymptotically exact 
in the low energy limit. 



A. Review on the Eliashberg theory 



We start from an effective U(l) gauge theory with one 
patch in one time and two space dimensions^ 2 - 



C = fl (i]d T - id x 
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where f a and a represent fermionic spinons and U(l) 
gauge fluctuations, respectively, emerging in the U(l) 
spin liquid stated, e is an internal gauge charge of 
a spinon and N is the spin degeneracy. r\ is an in- 
finitesimal coefficient to control artificial divergences in 
quantum corrections, which can be cured by self-energy 
corrections^, z is the dynamical exponent determining 
the dispersion relation of gauge fluctuations. It is given 
by z — 3 for several problems such as ferromagnetic or 
nematic quantum criticality including the present spin 
liquid problem 1 ^ while z — 2 in the spin density wave 
ordering 16 . Both the Fermi velocity vp and the curva- 
ture l/m are set to one. 

It was shown that either scattering with small momen- 
tum transfer or that with 2fcp (twice of the Fermi mo- 
mentum) is relevant in the Fermi surface problem 2 -^. In 
the one patch approximation^ only forward scattering, 
identified with 54 in the g-ology of the one dimensional 
problem 17 , is taken into account while another forward 
scattering (32) and back scattering (gi) are neglected. 
Such scattering channels can be introduced in the two 
patch approximation^. 

The previous large- N analysis without vertex correc- 
tions coincides with the Eliashberg approximation-ii, in- 
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troducing only self-energy corrections 
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where the spinon Green's function and gauge propagator 
are given by 

G a (k ,k) 



D(q ,q) = 



kx 
1 
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S(fco)' 
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respectively. The main point is that dynamics of gauge 
fluctuations is given by the Landau damping term with 
the damping coefficient 7^ proportional to kp , resulting 
from Fermi-surface fluctuations, while spinon dynamics 
shows non-Fermi liquid dependence in frequency of its 
self-energy, given by 2/z with a constant A;,. 

The problem to address in this paper is whether the 
anomalous exponent 2/z will be modified or not when 
vertex corrections are taken into account in a non- 
perturbative way, i.e., up to an infinite order. It turns 
out that the gauge dynamics cannot be modified from the 
Landau damping dynamics in the one patch approxima- 
tion while the fermion dynamics is expected to have some 
corrections^. Even in the two patch approximation, the 
gauge dynamics is still unchanged while fermions were 
shown to have 1 /N modification for the frequency expo- 
nent in the perturbative approach based on the Eliash- 
berg solution^. 
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FIG. 1: Fermion self-energy £(&) and boson self-energy 11(g), 
where the thick line represents the fermion Green's function 
and the wavy line does the gauge propagator. The shaded 
region can be any renormalized vertex. 
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FIG. 2: (Color online) The ladder vertex correction turns out 
to be irrelevant in the Eliashberg solution. 



B. Self-consistent ladder approximation 

We introduce vertex corrections in a completely non- 
perturbative way based on the ladder approximation. 
Then, we obtain full self-consistent equations 



n(go,q) = 




G a {k +q ,k + q)G a (k , k), 

s ( fc o) = ~^ J J ^4^A(k a + q a ,k + q:k ,k) 

G tT (ko + q ,k + q)D(q ,q), (4) 

where A(fco + qo, k + q; ko, k) is the vertex function given 
by 

1 1 e " f dl ° f rf2/ 

A(fc + ,o,fc + g;fc ,fc) = l-- j / ^ J 

A(/c + q - lo, k + q - I; k Q - l 0l k - l)D(l , 1) 
G a (k + qo - l ,k + q - l)G a (k Q -l ,k- I) (5) 



in the ladder approximation. Fig. 1 represents Eq. (4) 
and Fig. 2 displays Eq. (5). 

In order to solve these three coupled integral equations, 
we consider the Ward identityii 

iriq 1{k + q ,k + q: k , k) + q x A(k Q + q ,k + q; fc , k) 
+q y F(k Q + q ,k + q; k , k) 

= G- 1 (k + q ,k + q)-G- 1 (k 0: k) 7 (6) 

where ~I(fco + qo, k + q; ko, k) is the scalar vertex while 
A(fc + qo,k + q; k , k) and F(k + qo, k + q; k , k) are 
vector vertices. There is a special relation between the 
scalar vertex and vector one in one dimension due to the 
kinematic constraint, that is, the vector vertex is pro- 
portional to the scalar vertex with the Fermi velocity^. 
This is the reason why the one dimensional problem is 
exactly solvable. On the other hand, such a relation does 
not exist above one dimension, thus it is necessary to 
propose an ansatz for such a relation. 
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We suggests the following relation 
~l(k + q ,k + q;k ,k) -t A(k + q , k + q; k , k), 
F(k + qo,k + q;k a ,k) -> (2k y + q y )A(k a + q , k + q; k , k), 

(7) 

where the first ansatz is the application of the one dimen- 
sional result while the second one is our main assumption. 
As a result, the vertex function is 

: l (k + q ,k + q)-G- 1 (k ,k) 



The first part denoted by Si(fco) turns out to vanish 
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On the other hand, the second part recovers exactly the 
same expression as Eq. (2) in the low energy limit 
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where g~ 1 (ko,k) = iijko + k x + k 2 is a non-interacting 
Green's function. We note that this expression recovers 
both non-interacting and one dimensional cases. In ap- 
pendix A we show that this ansatz is self-consistent for 
the vertex equation [Eq. (5)] in the low energy limit. 

It is educational to check that if we apply the well- 
known one dimensional ansatz for the vertex function, 
i.e., neglecting F(k + q ,k + q; k , k), given by 

G- l {k + q ,k + q) 



where A is a modified constant, 
of this result in appendix B. 



We show the derivation 



Discussion 



A(fc + qo,k + q; k ,k) 

irjlo + q x 

we find that the Landau damping dynamics for gauge 
fluctuations is not reproduced as follows 
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This implies that the one dimensional Ward identity can- 
not be applied to higher dimensional cases. 

Inserting the vertex function Eq. (8) into the equation 
for the polarization function in Eq. (4), we find that the 
Landau damping dynamics does not change as follows 
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where 7 is a modified damping coefficient. 

Inserting both the vertex function [Eq. (8)] and boson 
self-energy [Eq. (9)] into the equation for the fermion 
self-energy in Eq. (4), we obtain two sectors 
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An essential point is that the ladder-vertex equation 
[Eq. (5)] is solved, resorting to the Ward identity, where 
the ansatz for the relationship between the scalar and 
G a (fcp,fc)v ector vertices [Eq. (7)] is introduced to result in the rela- 
tionship between the vertex function and fermion Green's 
function [Eq. (8)]. Justification for Eq. (8) lies in the 
fact that the boson self-energy should be given by the 
Landau damping solution. In the one patch formulation 
higher order quantum corrections are shown to vanish 
identically because all poles in the integral expression 
are in the same half plane, implying that the Landau 
damping solution is exact^. Considering the structure 
of the boson self-energy with the ladder- vertex correction 
[Eq. (9)], Eq. (8) seems to be generic. In addition, Eq. 
(8) recovers not only the non-interacting case but also 
the one dimensional physics completely. The relation- 
ship between Eq. (7) and Eq. (8) is unique as far as the 
y-current vertex is linearly related with the x-current or 
scalar vertex. 

Can we use the Ward identity in this approximation 
scheme? Usually speaking, the Ward identity is on the 
relationship between full vertex corrections and corre- 
sponding Green's functions. Actually, the Ward iden- 
tity and the special relation between the vector and 
scalar vertices in one dimension are satisfied for the 
fully renormalized vertex. In fact, the Ward identity 
should be always satisfied in any approximation scheme 
because it guarantees conservation of the system. In 
this sense the Ward identity may be regarded as an- 
other phrase of the conserving approximation. Math- 
ematically speaking, the conserving approximation can 
be derived from the Luttinger-Ward functional approach, 
where fully self-consistent sets of equations derived from 
the Luttinger-Ward functional respect the Ward identity 
automaticall y 18 ' 19 . Of course, this is not exact. A good 
example can be found in the impurity problem, called the 
conserving self-consistent t- matrix approximation 20 . 

Another important assumption is that singular depen- 
dence of the fermion self-energy occurs from frequency 
(10) instead of momentum. Though this is a common result 
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within the Eliashberg approximation^, there is no rea- 
son a priori, for this assumption to remain valid as soon 
as ladder-type vertex corrections are included. In par- 
ticular, the self-consistent calculations performed in the 
appendices are carried out within this assumption, which 
greatly simplifies the computations. It was argued that 
the fermion self-energy has the same frequency depen- 
dence as the Eliashberg solution and there is no singular 
momentum dependence in the perturbative evaluation of 
the one patch formulation up to an infinite order based 
on the Eliashberg solutio n 12 ' 14 , implying no anomalous 
exponent for the fermion Green's function, although one 
can not remove the possibility that the summation for 
coefficients from higher order quantum corrections may 
be singular. In addition, ladder-type vertex corrections 
turn out not to change the Elaishberg solution in the 
perturbative calculation of the two patch formulation up 
to the lowest order—. In this respect our result is not 
surprising but expected from the perturbative analysis 
in both one-patch 1 - 2 ' and two-patcb~ formulations. 

However, special types of quantum corrections involved 
with 2kp momentum transfer were shown to cause the 
singular momentum dependence for the fermion self- 
energy in the two patch formulation, giving rise to an 
anomalous exponent for the fermion Green's function 1 ^. 
Unfortunately, these quantum processes are not in- 
troduced in the ladder approximation, given by the 
Aslamasov-Larkin diagrams. This leads us to consider 
the Aslamasov-Larkin vertex up to an infinite order, 
shown in Fig. 3. Frankly speaking, this consideration 
is not completely new, already investigated in the con- 
text of the superconducting instability although it cor- 
responds to the particle-particle channel 2 ^. An interest- 
ing point is that such vertices are singulary enhanced 
in the U(l) spin liquid state, causing anomalous criti- 
cal exponents according to the perturbative evaluation. 
There is another Aslamasov-Larkin vertex correction in 
the particle-particle channel associated with supercon- 
ductivity, competing with the 2kp particle-hole instabil- 
ity. One problem in this consideration is to construct the 
self-consistent conserving approximation, not addressed 
clearly as far as we know. When this self-consistent con- 
serving framework is settled, we can check whether the 
vertex function [Eq. (8)] from the Ward identity allows 
the self-consistent solution with new critical exponents 
or not. If it works, we have a powerful framework. 

This discussion reminds us of the conserving self- 
consistent t-matrix approximation (CTMA) for the single 
impurity problem 2 ^, introduced to overcome the failure 
of the non-crossing approximation (NCA) in the exactly 
screened case, i.e., the absence of the strong coupling 
fixed point below the Kondo temperature. Here, the 
NCA is analogous to the Eliashberg approximation while 
the CTMA is parallel with the self-consistent 2kp treat- 
ment. Although the CTMA does not change the critical 
exponents of the NCA in the over-screened case, it cures 
several problems associated with thermodynamics^. In 
the present situation the self-consistent 2kp treatment 




FIG. 3: (Color online) The Aslamasov-Larkin vertex correc- 
tion in the particle-hole channel, where 2&f momentum trans- 
fer gives rise to an anomalous critical exponent for the fermion 
propagator. This situation is quite analogous to the super- 
conducting instability. Actually, there is another Aslamasov- 
Larkin vertex correction in the particle-particle channel, com- 
peting with the 2kF particle-hole instability. Inserting these 
Aslamasov-Larkin vertex corrections into the shaded regions 
of Fig. 1, one can construct self-consistent equations for 
fermion and boson self-energies. 



may change the critical exponents of the Eliashberg ap- 
proximation. It will be really interesting to investigate 
this type of diagrams near future. 



III. SUMMARY 

In this study we solved three coupled integral equa- 
tions for boson and fermion self-energies with vertex cor- 
rections. The key point is to make an ansatz for the 
vertex function in terms of the fermion Green's function 
based on the Ward identity. Resorting to this ansatz, 
we find a set of full self-consistent solutions in the ladder 
approximation for the vertex correction, where both bo- 
son and fermion self-energy corrections do not have any 
modification, compared with the Eliashberg approxima- 
tion. This implies that if the anomalous exponent arises 
in the fermion self-energy, it may result from the class 
of Aslamasov-Larkin diagrams, not taken into account in 
our study. This remains as an important future work. 
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vating us to investigate the present problem. K.-S. Kim 
was supported by the National Research Foundation of 
Korea (NRF) grant funded by the Korea government 
(MEST) (No. 2010-0074542). 



Appendix A: Self-consistency between the Ward 
identity and the ladder approximation for the vertex 
function in the low energy limit 

In appendix A we show that the ladder vertex in Eq. 
(5) satisfies the Ward identity with the ansatz Eq. (8). 
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In other words, the ansatz Eq. (8) is shown to be self- 
consistent at least in the ladder approximation. 

Inserting Eq. (8) into the right hand side (R.H.S) in 
Eq. (5), we obtain 



G a (k -l ,k-l) - G CT (fc + q -lo,k + q-l) 



9a 1 (k + q -lo,k + q-l)- ga 1 (k a -l ,k- I) 
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Performing the l x integration, we obtain 



Expanding the first sector in l y , we obtain 
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Integrating over the frequency Iq, we obtain 
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The first In contribution can be approximated as follows 
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One will realize that the first term is associated with the 
self-energy in the Eliashberg approximation. Evaluating 
the B term in the same way as A and gathering both A 
and B, we reach the final expression 



A(fc + q , k + q; fc , k) - 1 = A + B 

S(fc + gp) - S(fcp) 

9a 1 (ka + q ,k + q) - ga 1 (k Q ,k) 
+J r (k + q ,k + q;k Q ,k), (A6) 
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It is not difficult to observe that J"(fco + qo, k + q; ko, k) 
is irrelevant in the low energy limit due to the frequency 
and momentum dependence in the numerator, giving rise 
to higher order corrections to the fermion self-energy. We 
conclude that Eq. (8) is asymptotically correct in the low 
energy limit. 
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Appendix B: Evaluation of £2(^0) The vertex correction does not change the scaling for 

frequency in the ladder approximation. 

In appendix B we evaluate £2(^0)- Performing mo- 
mentum and frequency integrals, we obtain 
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dq i-qko + k x + k 2 - S(fc ) 
2n ir)k + k x + k 2 - S(fc + g ) 
sgn (go )- i Sgn( fco + go ) ^ A sgn(fco)|fco|2/ ,_ (m) 
% 
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